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High-order harmonic generation (HHG) in simultaneous intense near-infrared (nir) laser light 
and brilliant x rays above an inner-shell absorption edge is examined. A tightly bound inner-shell 
electron is transferred into the continuum. Then, nir light takes over and drives the liberated 
electron through the continuum until it eventually returns to the cation leading in some cases to 
recombination and emission of a high-harmonic photon that is upshifted by the x-ray photon energy. 
We develop a theory of this scenario and apply it to Is electrons of neon atoms. The boosted high 
harmonic light is used to generate a single attosecond pulse in the kiloelectronvolt regime. Prospects 
for nonlinear x-ray physics and HHG-based spectroscopy involving core orbitals are discussed. 

PACS numbers: 41.60.Cr, 32.80.Aa, 32.30.Rj, 42.65.Ky 



High-order harmonic generation (HHG) by atoms in 
intense near-infrared (nir) laser fields is a fascinating 
phenomenon and a versatile tool; it has spawned the field 
of attoscience, is used for spectroscopy, and serves as a 
light source in many optical laboratories [l|-|3|. Within 
the single-active electron (SAE) approximation [^-|6[ and 
the three-step model of HHG [Tl, |8| , the NiR laser tun- 
nel ionizes a valence electron and accelerates it in the 
continuum. When the nir laser field changes direc- 
tion, the liberated electron is driven back to rescatter 
with the parent ion. This may cause the electron to rc- 
combine with the ion whereby the excess energy due to 
the atomic potential and due to the energy gained from 
the NIR laser field is released in terms of a high har- 
monic (HH) of the NIR frequency. Frequently, this mind- 
set is also applied to two-color HHG where a NIR laser 
is combined with vuv/xuv light [1, [13] which thereby 
assists in the ionization process leading to an overall in- 
creased yield [9|, [llj, [l^ . This principle is evolved fur- 
ther by using attosecond xuv pulses to boost the HHG 
process which increases the yield for a certain frequency 
range by enhancing the contribution from specific quan- 
tum orbits ISl - llSl . However, there are only few excep- 
tions, e.g., Refs. ll9l422l . in which many-electron effects are 
treated for two-color HHG. Fleischer [20| includes implic- 
itly other electrons by using a frequency-dependent po- 
larizability for the atoms. Then the xuv light is found 
to cause new plateaus to emerge at higher energies, how- 
ever, with a much lower HH yield. Explicit two-electron 
effects in two-color HHG are considered in Refs. 21 and 



In this letter, we consider even higher photon energies 
where core electrons directly couple to the continuum 
using x rays from a free electron laser (FEL) such as 
the Linac Coherent Light Source (LCLS) [23J, |2J]. Our 
SAE scheme is depicted in Fig. [T] and proceeds in allusion 
to the three-step model of HHG [7|, ISJ as follows: (a) a 
core electron is ionized by one-x-ray-photon absorption; 
(b) the liberated electron propagates freely; (c) in some 
cases the electron is driven back to the ion and recom- 
bines with it emitting HH radiation. Although Fig. [1^ 
suggest that the electron is born at a specific phase of 
the NIR light, this is not the case; instead, electrons are 
ejected during the entire x-ray pulse. This lifts restric- 
tions of the width of the HHG plateau similarly to stud- 
ies of HHG with valence- ionization by xuv light |16l4l8| . 
We exemplify our method for Is electrons of neon with 
an ionization potential (IP) of h = 870.2 eV ^. This 
yields boosted HH radiation at close to kiloelectronvolt 
photon energies from which we isolate a single attosecond 
pulse. All equations are formulated in atomic units. 

Our scheme offers completely new prospects for HHG 
by involving core electrons. It thus goes beyond just 
an extension of the HHG cutoff into the kiloelectronvolt 
regime for which one may use conventional concepts that 



221 : there, the xuv photon energy is tuned to the core- 
valence resonance in the transient cation that is produced 
in the course of the HHG process by tunnel ionization. 
This leads to a second high- yield plateau that is shifted 
to higher energies by the xuv photon energy. 





(a) 



FIG. 1. (Color online) Schematic of the three-step model for 
the HHG process modified by x-ray induced ionization of a 
core electron. See text for details. 



are laboratory size and valence-electron-based. Here it 
has been shown recently using midinfrared lasers [26[ at 
high intensities of ~10^^ W/cm^ that good HH yields can 
be obtained even in the presence of enhanced ionization 
by a judicious choice of the gas pressure |27H29| . 

In our quantum theory of x-ray ionization-based HHG, 
we make the SAE \^]^ . The ground state | ) is repre- 
sented by the electron in the core orbital of the atomic 
ground state. Continuum states are described by plane 
waves I fc ) for fc G R'^ where we neglect the impact of 
the atomic potential which is called the strong-field ap- 
proximation. There arc two processes which destroy the 
system and inhibit HHG. First, even in field-free condi- 
tions, core holes decay with a corresponding width Tc- 
Second, for two-color light additional valence-shell pho- 
toionization occurs for neutral and core-ionized neon in- 
duced by both, the nir laser and the x rays [30[. The 
total destruction rates are expressed as ro(i) and T^{t) 
for the neutral and the core-ionized states, respectively. 
They depend on time because the photoionization rate 
depends on the envelope of the nir and x-ray pulse. Here 
"A:" is only a label as ^^{t) depends only negligibly on k 
for fast continuum electrons. 

Based on the above model of an atom in two-color light, 
we derive equations of motion (EOMs) for its time evo- 
lution using the following ansatz for the SAE wavepacket 
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where /p is the IP of the core electron. The Hamiltonian 
of the model is H = Ha + Hi^ + Hx', it consists of the 
atomic electronic structure Hamiltonian Ha and the cou- 
pling of the SAE to the nir laser Hi^ and the x rays Hx in 
electric dipole approximation in length form [31[ . Insert- 



ing Eq. ID) into the time dependent Schrodinger equation 
and projecting onto (0| yields the EOM for the ground- 
state amplitude 



iait) = --roit) ait) + Exit) I bik,t){0\ex-r\k) d'^k 



Projecting onto ( fc | for all k G 
the continuum amplitude 
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results in EOMs for 
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where r is the electric dipole and ex is the linear po- 
larization vector of the x rays. Inner-shell electrons 
are tightly bound such that the nir laser hardly affects 
them. Hence, in the transition matrix elements only 
the x-ray term is relevant, i.e., {k\Hx |0). Conversely, 
only the nir laser impacts continuum electrons noticably, 
i.e., (fc|i?L|fc')=iVfc<53(fc_fc'). 



To solve the coupled system of first-order partial differ- 
ential equations ([2]) and (jSj, we realize that we may ne- 
glect the second term on the right-hand side of Eq. ^ [5[ . 
This approximation decouples Eq. ^ from Eq. ([3|) which 
can now be integrated directly assuming that there is no 
light for t < 0. Then, the solution is unity for i < and 
a(t) = Q-2Fo{t) for t > where the temporal destruc- 

t 
tion exponent is defined by F iit) = jViit') dt' for i G 



{0}U{A:|VfcGM^}. With the closed-form solution of 
Eq. (121), we can now integrate Eq. ^ formally exactly 
by introducing the canonical momentum p ~ k — Ai^ it) 
with the vector potential of the nir laser ^l(0 I^- 
The saddle point approximation is used to simplify the 
triple integration over p in the calculation of the electric 
dipole transition matrix element Dit) = J a* it) (0| eb ■ 

f\k) bik,t) d^k ~ £)'(t)e~'"''* aloiig the direction eb 
which determines the HH emission 0, 0, S, [iH, [H, ^ 
with the slowly varying dipole moment 
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Wc introduced the excursion time t = t — t'\ at the sta- 
tionary point (saddle point), the momentum is Psti^i t) = 

1 * - 

/ Ai^it') dt' and the quasiclassical action is 
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Further, we employed the rotating- wave approxima- 
tion [31| such that in Eq. ^ only the positive frequency 
components of the x-ray electric field E^it) are taken 
leading to the term involving "— wx" in Eq. ([5]) via the 

decomposition E^it) ~ — Exit) e~'"^* with the complex 

field envelope Exit) = Eoxit)c-'^^'''-*^+^''-"\ the real 
field envelope Eoxit), the time-dependent phase (fixit), 
and the carrier to envelope phase (CEP) ifx.o |33l |. 

The HHG spectrum follows from Dit) by Fourier trans- 
formation. The transforms of -D(i) and D'it) are related 
by Diuj) = D'iuj—LUx), i.e., the entire harmonic spectrum 
is shifted by wx toward higher energies. The harmonic 
photon number spectrum (HPNS) [2l|, |22|, |32| of a single 
atom — i.e., the probability to find a photon with speci- 
fied energy — along the propagation axis of the NIR laser 
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FIG. 2. (Color) Two exemplary SASE x-ray FEL pulses 
(black and red) for an x-ray photon energy of uy^ — 569 ojl ~ 
881.8 eV. The envelope Eox{t) is specified in units of the 
peak electric field _Egx for the given peak x-ray intensity 
of 7.5 X 10^* W/cm^ and is represented by the solid lines 
whereas the phase ipx{t) is depicted as dashed lines with 
CEP (px,o ~ 0. The average FWHM pulse duration is 1.5 op- 
tical cycles with a cosine square pulse envelope [3J| . The aver- 
age pulse spectrum is Gaussian and has a FWHM bandwidth 
of 8eV which corresponds to a coherence time of 0.24 fs. 



with the density of free-photon states g{uj) ~ ,^^3 ^3 [3l|, 
the speed of light in vacuum c, and the sohd angle fl. 

We apply our theory to generate HH radiation from 
Is core electrons of a neon atom where the polariza- 
tion vectors of nir laser bl and x rays ex and eb are 
along the z axis. To calculate the Is orbital of neu- 
tral neon, we use the program of Herman and Skillman 
setting a — 1 [35l437f . For the experimental core- hole 



width, Fc = 0.27 eV, is used from Ref. |25|. Present-day 
FELs generate x rays on the basis of the SASE prin- 
ciple [2j, I2J] . We model SASE pulses with the partial 
coherence method J38l - l4d |. Two sample pulses are dis- 
played in Fig. [2] where modeling parameters are specified 
in the figure caption. The chosen wx = 881.8 eV is above 
the K edge but still close to it such that ionization by 
X rays is very efficient 
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The rate of destruction of 
neutral and core-ionized neon by the nir laser and the 
X rays is determined as follows. For photoionization by 
X rays, we use the total valence cross section (Ne 2s and 
Ne2p) at utx which is ctq = 2.3 x 10' 
tral atom and ctj: = 3.1 x 10" 
cation obtained with Refs. |42| and l43l For the NiR laser, 
destruction by tunnel ionization occurs with the instan- 
taneous rate Fo,L(i) that is determined by the ADK for- 
mula 13, S] and the valence IP of neon 21.5645 eV J46l |: 
at the chosen peak intensity [caption of Fig. |3] , the cycle- 
averaged destruction rate is 0.009 eV. The ionization rate 
for the cation is set to zero, i.e., rjr^(t) = OeV, be- 
cause the valence IP of core-ionized states is much larger 
than the valence IP of the neutral atom leading to a 

vanishingly small width due to the exponential depen- 

I — I PHi 
dence of the ADK rate on the IP [4J, |45| . Summing up 

all contributions yields Tf){t) = ao Jx{t) + To^hit) and 

Tfiit) = Fc -l-CTj: Jx(i) +Fj: ^(i) where Jx(i) is the instan- 
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FIG. 3. (Color online) Harmonic photon number spec- 
tra (HPNS) [Eq. de])] (a) for Ne 2p valence electrons and (b) for 
Ne Is core electrons using the SASE pulses of Fig. [21 The 
solid, black line is from the black SASE pulse and the dashed, 
red line from the red pulse. The x-ray intensity [caption of 
Fig. [2] was chosen such that the ionization rate is the same 
as the tunnel ionization rate by the nir laser. The nir laser 
vector potential A-L{t) has a cosine square pulse envelope [SJI 
with a FWHM duration of 1.5 optical cycles, CEP ^l,o = 7i"/2, 
and a peak intensity of 3 x 10^* W/cm^ for 800 nm central 
wavelength, i.e., the central photon energy is ojl = 1.55 eV. 



taneous x-ray flux. 

HHG spectra from Ne Is core electrons for the two 
SASE pulses of Fig.[2]are shown in Fig.[3]D together with a 
conventional HHG spectrum from Ne 2p valence electrons 
in Fig. [3^ to facilitate a comparison. In Fig.[3]D, we see an 
extension of the spectrum toward higher harmonic orders 
as well as a significant impact of the SASE pulse shape 
on the spectra as was also found for x-ray boosted HHG 
based on core excitations [2lJ, |22|. To gain deeper insight 
into the temporal evolution of the emission of HHG light 
from core electrons, we carry out a time-frequency anal- 
ysis of the dipole moment D{t) giving D{t,uj). This is 
done with a windowed Fourier transform using a Gaus- 
sian window with a variance of 0.1 fs. The result corre- 
sponding to the black SASE pulse of Fig. [3] is shown in 
Fig. 21 Comparison of this plot with the time-frequency 
analysis of the red SASE pulse of Fig. [2] and the inves- 
tigation of a difference plot reveals that the imprinting 
of the SASE pulse onto the HHG emission is noticeable 
but the general structure is determined by the nir laser. 
For the 1.5-cycle nir field, the intensity for each cycle is 
different, causing the highest harmonic orders to be pro- 
duced at the peak of the NiR laser pulse. Then, the single 
maximum of the nir laser vector potential produces the 
harmonic orders larger than 640 wl . By filtering out har- 
monic orders lower than 640 wl, we obtain the single at- 
tosccond pulses for the two SASE pulses that wc fit with 
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FIG. 4. (Color) Time-frequency analysis of core-electron 
HHG in neon for the black SASE pulses of Fig.[l 



a Gaussian profile of FWHM duration of about 160 as. 
The SASE pulse shape impacts only the peak intensity 
of the attosecond pulses. 

To evaluate the efficiency of x-ray boosted HHG with 
Ne Is core electrons, we compare it with the nir laser- 
only HHG with Ne2p valence electrons [Fig. [3]. Both 
X rays and nir laser give rise to the same peak ioniza- 
tion rate. Inspecting Fig. [31 we find that the efficiency of 
x-ray boosted HHG is reduced by a factor of about 10~^ 
compared with the case of HHG with only the nir laser. 
This can be explained by the following differences be- 
tween the two HHG processes. First, the initial ioniza- 
tion step is mediated by either tunnel ionization or a 
one-x-ray-photon absorption process. Tunnel ionization 
occurs predominantly along the linear NIR laser polar- 
ization axis leading to a continuum electron wavepacket 
elongated along this axis. In contrast, one-x-ray-photon 
ionization of a Ne Is electron produces a wavepacket with 
the figure-8-shape of a pj,-orbital. Hence only the part of 
the electron which is ejected along the nir laser polar- 
ization axis will have a significant probability to recom- 
bine with the ionic remnant and thus contribute to HHG. 
Second, core-hole recombination is less efficient than 
valence-hole recombination. Third, valence holes are sta- 
ble whereas core holes decay. Fourth, tunnel ionization 
releases electrons at rest but one- x-ray-photon ionization 
ejects electrons with a kinetic energy of wx — ^p |4ll |. 

Finally, we need to investigate phase matching of 
x-ray boosted HHG in a macroscopic medium, which is 
crucial for predicting the output of HHG, but becomes 
increasingly difficult with rising HH photon energies 27 - 
29l . |47| . From Eq. ([6|) , we estimate the emitted photon 

number per pulse by iVph — {rig, Vy Afl J ^^^^^ dw, 



waist Wo, the Rayleigh range L = ttwq/Al, and the 
NIR laser wavelength Al [11] ■ Next AD, « 2 7rAi?^ is 
the solid angle of the cone of HH emission that has an 
opening angle d along the propagation axis. We have 
the condition ^^ < \i^/{Lh) for constructive interfer- 
ence in terms of the highest-generated harmonic order h. 
When h = 660, wq = 30 Al, and n^ = 10^^ cm'^, 
then iVph is estimated by an order of magnitude to be 
about 700 which is a sufficiently high to make our boosted 
HHG light interesting for applications. In the attosec- 
ond pulses there are about 30 photons. Furthermore, 
A^ph can be increased by orders of magnitude using a 
higher x-ray intensity and — if one is only interested in a 
high HH yield — longer light pulses. 

In conclusion, using x rays to ionize core electrons as 
a first step in HHG, we predict the high harmonics be- 
ing boosted into the x-ray regime. By superimposing the 
harmonics close to the cutoff, we have obtained single at- 
tosecond pulses. Our scheme extends most HHG-based 
methods to inner-shell atomic and molecular orbitals us- 
ing suitably tuned x rays. 
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